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> ■ Abstract 

^: 

l/^ . We develop a microscopic analysis of superconducting and dissipative cur- 

^^ i rents in junctions between superconductors with d-wave symmetry of the 

,-_i I order parameter. We study the proximity effect in such superconductors and 

^ ' show that for certain crystal orientations the superconducting order parame- 

f— ^ ' ter can be essentially suppressed in the vicinity of a nontransparent specularly 

^—^ ■ reflecting boundary. This effect strongly influences the value and the angular 

I ! dependence of the dc Josephson current js- At T ~ Tc it leads to a crossover 

■^ I between js oa Tc — T and js oc {Tc — T)^ respectively for homogeneous and 

^^ ' nonhomogeneous distribution of the order parameter in the vicinity of a tun- 

C^ ' nel junction. We show that at low temperatures the current-phase relation 

H . jsiv) for SNS junctions and short weak links between d-wave superconduc- 

— L , tors is essentially nonharmonic and contains a discontinuity at c/? = 0. This 

Ch ', leads to further interesting features of such systems which can be used for 

j^ I novel pairing symmetry tests in high temperature superconductors (HTSC). 

kJ ' We also investigated the low temperature I-V curves of NS and SS tunnel 

• '-j ■ junctions and demonstrated that depending on the junction type and crystal 

i^j] , orientation these curves show zero-bias anomalies / oc V"^, I oc y^ln(l/l/) 

P^ \ and / oc V^ caused by the gapless behavior of the order parameter in d-wave 
superconductors. Many of our results agree well with recent experimental 
findings for HTSC compounds. 

PACS numbers: 74.50. -hr, 74.20.-z 
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I. INTRODUCTION 

In spite of enormous efforts made to understand the physical mechanisms of pairing in 
various high temperature superconductors (HTSC) the situation still remains unclear. A 
key role in understanding of this phenomenon belongs to the question about the symmetry 
of the order parameter. Quite early after the original Bednords and Miiller discovery |l[] the 
symmetry of the dx2-y2-tYpe was suggested for HTSC materials 0-^. Since then a plenty of 
experiments have been designed to probe the symmetry of the order parameter in HTSC (see 
e.g. 10 for a review). Although many experimental results are consistent with the picture 
of d-wave pairing (e.g. the temperature dependence of the penetration depth 0], NMR and 
NQR studies |§] etc.) they still do not allow to rule out other possibilities, like anisotropic 
s-wave pairing. Moreover the results of some other experiments (see e.g. 0) may indicate 
s-wave rather than d-wave symmetry of the order parameter in HTSC. Therefore it is quite 
likely that only a set of different and independent experimental tests would allow to make 
an unambiquous conclusion about the order parameter symmetry in HTSC compounds. 

An important information about the symmetry of superconducting pairing can be ob- 
tained from the measurements of both the dc Josephson effect and the quasiparticle current 
in tunnel junctions between two HTSC. The dc Josephson effect in unconventional supercon- 
ductors has been discussed by Geshkenbein, Larkin and Barone |]TU[ and by Sigrist and Rice 



TT| who demonstated that the d-wave symmetry of the order parameter may lead to the 
sign inversion of the Josephson critical current for certain crystal orientations. Under these 
conditions the tunnel junction becomes the so-called vr-junction [|12|. Being closed by a vr- 
junction a SQUID loop with a not very small inductance develops a spontaneous circulating 
current [|l^]. As a result the magnetic flux equal to a half of the flux quantum occurs inside 
a ring and can be easily measured. Measurements of that kind have been carried out for 
HTSC samples [|13]-[T^ and indeed demonstrated the results fully consistent with the above 
picture. These results in combination with the paramagnetic behavior of granular HTSC 
compounds |T6[ and its theoretical interpretation [pA] , p!7| , p!8| serve as a serious argument in 



favour of d-wave pairing symmetry in HTSC. 

In contrast to the Josephson effect which is sensitive to the order parameter phase dif- 
ference across the junction low temperature measurements of the quasiparticle current in 
tunnel junctions provide information about the quasiparticle density of states in supercon- 
ducting banks and allow to distinguish gapless superconductivity from that with a finite 
gap. The results of numerous experiments vary from a nearly BCS-like to a clear gapless 
behavior for different HTSC materials (see e.g. and references therein). The results of 
recent tunneling experiments with Bi2Sr2CaCu208 samples 0,|3l indicate a gapless behav- 
ior of the I-V curve at low voltages and temperatures (e.g. the dependence J oc \^^ at low 
V was reported in [p!9| ). 

A growing number of experimental data makes it necessary to develop a detailed analysis 
and specify theoretical predictions concerning both dc Josephson effect and quasiparticle 
tunneling between d-wave superconductors. Several phenomenological calculations assuming 
d-wave pairing have been already done (see below) attempting partial understanding of some 
experimental results. Nevertheless a number of important questions still has to be adressed 
in this context. E.g. a uniform distribution of the superconducting order parameter on 
both sides of a tunnel barrier was assumed in many calculations. Being obviously correct 



for isotropic s-wave superconductors this assumption may fail for d-wave ones depending 
on their orientation relative to the tunnel barrier plane. Below we will show that the 
spatial dependence of the order parameter due to the proximity effect becomes particularly 
important close to the superconducting critical temperature T^ having a strong impact on 
the dc Josephson effect in d-wave superconductors. 

Another problem appears if one applies the tunneling Hamiltonian method to investiga- 
tion of the charge transport through tunnel junctions in the d-wave case. The momentum 
dependence of the matrix elements describing tunneling between superconductors becomes 
particularly important in this case. It is easy to show that the choice of the tunneling ma- 
trix elements as being independent on the momentum direction (standard for the s-wave 
case) leads to confusing results for d-wave superconductors. Furthermore, an unambiguous 
choice of this dependence cannot be done within this method. This emphasizes the necessity 
to provide a microscopic description of tunneling between such superconductors based on 
matching of the electron propagators at the tunnel barrier. This approach leaves no space 
for ambiquity and, on top of that, it is not confined to the case of low transparency barriers 
but allows to study other types of weak links with highly transparent interfaces. 

In this paper we will provide an extensive microscopic study of the charge transport in 
various types of junctions between d-wave superconductors with BCS-like behavior of the 
density of states. The paper is organized as follows. In Section 2 we develop a detailed study 
of the proximity effect for d-wave superconductor-insulator and d-wave superconductor- 
normal metal structures. We investigate the spatial dependence of the superconducting 
gap function for various crystal orientations and temperatures and show that for particular 
orientations the gap at the superconductor-insulator interface can be completely suppressed. 
The dc Josephson current through tunnel junctions, SNS junctions and short weak links 
between d-wave superconductors is examined in Section 3. Our analysis allows to discover 
several new qualitative features of the Josephson effect in such systems which can be used 
for further experimental tests of the pairing symmetry in HTSC In Section 4 we investigate 
the I-V curves for superconductor-superconductor (SS) and normal metal-superconductor 
(NS) tunnel junctions. For most of crystal orientations we found gapless nonOhmic behavior 
in the limit of small voltages and T = 0. Discussion of our results is presented in Section 5. 

II. PROXIMITY EFFECT IN D-WAVE SUPERCONDUCTORS 

The order parameter in bulk superconductors with unconventional pairing depends on 
the direction of the Fermi momentum pp [pl|p2| . Beyond that close to the edges of a 
superconducting piece of metal the order parameter acquires a spatial dependence due to 
the proximity effect. In this section we present a detailed investigation of this effect for 
superconductors with d-wave symmetry of the order parameter. We show that the spatial 
dependence of the order parameter in the vicinity of a low transparency insulating barrier 
may essentially depend on the crystal orientation relative to the barrier plane. Similar 
- although quantitatively different - results hold provided a superconductor is in a good 
electric contact with a normal metal. 

In order to describe the proximity effect in d-wave superconductors we make use of the 
Eilenberger equations for the quasiclassical Green functions 123]. In the case of supercon- 



ductors with singlet pairing these equations read [p3|,p4 



r {2ujm + vpVRJfiP, R, ujm) - 2A(p, R)g{p, R, ujm) = 0, 

I (2cu„ - VFVR)f+{p, R, Um) - 2A*(p, R)g{p, R, Um) = 0, 

[ VF^RQiP, R, cu„) + A(p, R)f+{p, R, Um) - A*(p, R)f{p, R, Ur, 



(1) 



0. 



Here Um = (2m + l)7iT is the Matsubara frequency, p = Pf/\pf\, Vp{p) = pp/m is the 
Fermi velocity, A(p, R) is the order parameter or the gap function. Anomalous and normal 
Green functions f{p,R,uJm) = f^*{~P, R,^m) and g{p,R,uJm) = g*{—p,R,uJm) obey the 
normalization condition 



5f^(p, R, Um) + /(P, R, UJm)f^{P, R, UJm) = 1- 



(2) 



The order parameter A is linked to the anomalous Green function by means of the standard 
selfconsistency equation 
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where V{p,p') is the anisotropic pairing potential and the integration is carried out over 
the Fermi surface. The quasiclassical equations (|I]) are valid at the scale much larger than 
the interatomic distance ~ 1/pf and do not keep track of rapid changes of the system 
parameters very close to the metal-metal or metal-insulator boundaries. In order to take 
these boundary effects into account the system of equations (|l|)-(||) should be supplemented 
by the boundary conditions [P5|-p7[| matching the quasiclassical electron propagators g and 
/ on both sides of the boundary. These boundary conditions may essentially depend on 
the quality of the interface. In the case of a nonmagnetic specularly reflecting boundary 
between two metals these conditions read [p^ : 



d-{p_) = d+(p+). 
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Here [a, b] denotes the commutator of matrices a and b, R{p) is the reflectivity coefficient 
and the index +(— ) labels the electron momentum in the right (left) halfspace with respect 
to the boundary plane. The 2x2 matrices d and s are defined by the equations d{p) = g{p) — 
g{jp),s{p) = g{p) + g{p), where p{p) denotes the incident (refiected) electron momentum 
and 
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Provided the transparency of the tunnel barrier is equal to zero D = 1 — R 
(I) yield m 



(5) 



the equations 



9±ip)=9±{p), f±{p) = f±{p), mp)=f±ip), 



(6) 



where the Green functions are taken at the metal-insulator boundary. In the opposite 
limiting case of a transparent boundary between two metals D = 1 the equations @ reduce 
to a simple continuity conditions for the Green functions (7+ = g^ and /+ = /_ at the 
boundary plane. 

For the sake of definiteness let us assume that a d-wave superconductor occupies a half 
space a; > 0. Provided there is no current flow in this superconductor one can choose the 
gap function A to be real there and define /i = (/ + /'*')/2, /2 = (/ — /^)/2- Then with 
the aid of (0) we find 

?? /\ ( T} T] 1) 

f\{p, X, UJm) - -Ardlfiip, X, UJ^) - , , ' {I + -A^idJiip, X, iVm)? " /'(p, X, uj}^/^ = 0, 



^OOl '^■'''-' ' ■"' Urr, ^ AU. 
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(7) 

Let us first consider the case of an impenetrable boundary situated at the plane x = 0. 
Then the spatial dependence of the gap function A(x) is defined by the combination of 
the equations (P, (^ with the boundary conditions @ at x = and /i(x — > oo) = f^o, 
A{x -^ oo) = Aoo, where /oo and Aqo are the equilibrium values for the anomalous Green 
function and the order parameter in the bulk superconductor. Provided the order parameter 
obeys the condition 

A{p,R) = AiP,R) (8) 

the solution of the above equations does not depend on x and thus the functions /i and 
A coincide with their equilibrium values far from the boundary. For superconductors with 
d-wave symmetry of the order parameter A{pf) of the {pl^ ~P^o)"^yP^ ^^^ condition (||) is 
satisfied if one of the principal crystal axes Xq, yo or ^o is perpendicular to the boundary 
plane. 

For other crystal orientations the order parameter A{p,x) turns out to be spacially 
inhomogeneous. To proceed further let us assume that one can express the value A in 
the form A{p,x) = ip{p)r]{x) [^. At temperatures close to T^ and distances larger than 
the correlation length C,o from the boundary the function r]{x) obeys the Ginzburg-Landau 
equations which have a well known solution 

Vgl{x) = 7]^tanh[{x + (3)/V2C{T)], (9) 

r/oo is the equilibrium value of ri{x) far from the boundary and ^(T) is the temperature 
dependent superconducting coherence length. In the case of uniaxial symmetry we have 
^(T) = (^ji (T) cos^ a+^^(T) sin^ aY^'^, where ^||(T) and ^±{T) are the values of the coherence 
length respectively in the basal plane and the transversal direction, a is the angle between 
the vector normal to the boundary and the basal plane. The value of (3 is defined by the 
boundary condition 

qr]'{0) = 7/(0) 

and the parameter q has to be derived from the microscopic theory (Eqs. (|^), (^) and (|^). 
In the vicinity of the critical temperature T r^ Tc one can linearize the equation (^ 
neglecting higher powers of /i. Then combining (|^) and (|^) one gets 



fi{p,x,ujm) = 1 — r / {exp(-| — -{x-x')\)A{p,x') +exp(-| — -\{x + x'))A{p, x'))}dx' . 

\Vx\ Jo Vx V^ 

(10) 
Substituting ([10|) into (^ and setting A(p, x) = ip{p)rj{x) we arrive at the integral equation 

+ ^exp(-|^|(x + x'))W. (11) 

\Vx\ V^ 

The effective couphng constant for an anisotropic superconductor A is defined by the equation 

mp) = -/^(P'P')^(p')(^^' (12) 

which yields TfTcXJ^ml^ml^^ = 1- The equation (|n|) describes the behavior of ri{x) at 
distances x ^ ^(T) from the boundary. This equation coincides with that derived in ||3^ 
within the framework of a different technique for the case of a small gap anisotropy. 

A trivial combination of the above equations also allows to evaluate the Green functions 
at T close to Tc. E.g. the functions /i and /2 at the superconductor-insulator boundary 
X = and T —>■ Tc read 

r .. X V^(P)+^(P) r ( \'^^^\ \ I \A 
Ii{P,'^m) = , — , / exp(- \x)r]{x)dx, 

j2[P,^m) = sgnum exp(-| \x)r]{x)dx (13) 

Vx Jo Vx 

Ipip) - Ipip) 

= //Ml //-N Sgn(^x^m)/l(p,^m)- 
V^ P + V^ P 



The exact solution of (11) can be easily found for the case '?/'(p) = —ijj{p). Then we have 

ri{x) = Cx, 

i.e. for this particular crystal orientation the gap function vanishes at the boundary x = 
and for any x > it is described by the function 

ri{x) = ri^tanh[x/V2^{T)]. 

Combining the equation ip{p) = —iplp) with the symmetry condition for the order parameter 
of the {pI^ — Py^ytype we come to the conclusion that the gap function vanishes at the 
superconductor-insulator interface provided the principal crystal axes xq and yo constitutes 
the angle it/A with this interface. Note that this conclusion is essentially based on the 
symmetry arguments and thus remains valid at any temperature below T^. Indeed, for 
a pairing potential with the symmetry property V{p,p') = —V{p,p') and the boundary 
condition f{p) = f\p) we obtain 77(0) = from the selfconcistency equation (^ at any T. 
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At distances x ^ ^o from the boundary the integral equation ([TlD has a general solution 

r]{x) =C{x + q). (14) 

The parameter q can be easily found from a simple variational procedure [^ which yields 

' - U36C(3)r l>o^^^^^ + ^^^^J "^-^ ^ + 4vr3T /,^,,[^(p) - ^^(p)]^..^^^ 



^^(p)t;id^5)-\ (15) 

where integration over the Fermi surface is confined to its part with f^; > 0. 

Note that with the aid of general symmetry arguments one can unambiguously fix only 
those crystal orientations for which the parameter q is equal to zero and infinity. The detailed 
form of g as a function of a crystal orientation relative to the boundary plane for a given 
pairing symmetry type essentially depends on the particular choice of the basis function 
ip{p). E.g. for the pairing symmetry of the type {p^^ — Py^) for the sake of definiteness one 
can choose the simplest basis function '?/'(p) = Aq{pI^ ~'Pyo) ^"^^ after an explicit integration 



in (|T5|) get 

4C(3)t;^ r (cos^ 00 + 3 sin^ 00 cos^ 



^o)^ 



^ 57r3Te(l + 2 sin^ Oq) I sin^ 0o(l - sin^ Oq cos2 20o) 

+ ^^Q)'(4cos2 0o + 19sin^0oCos220o)}. (16) 

Here Oq, (po are the polar and the azimuthal angles of the normal n to the boundary (see 
Fig.l). The function ^(^o, v^o) is presented in Fig. 2. We see that the parameter q varies from 
q = {or T] = 0) for Oq = 7r/2, 0o = tt/A {ip{p) = —ip{p)) to g = cxd (i.e. ri'{0) = 0) if one of 
the main crystal axes Xq, Hq or zq has the angle 7r/2 with the interface plane. 

Our analysis can be easily modified to take into account the effect of nonmagnetic im- 
purities. In the presence of such impurities close to Tc the equation (pJ^ remains valid 
if one substitutes Um — ^ ^m = ^m + sgnc<Jm/2rjmp, where Timp is the average scatter- 
ing time. Accordingly the modified critical temperature T'^ is defined by the equation 
ttT^'A X]-m Pm|~^ = 1- With the aid of (|ll]) it is easy to check that scattering on nonmagnetic 
impurities does not lead to qualitative changes in the behavior of the order parameter in the 
vicinity of a superconductor-insulator interface. E.g. the homogeneous solution ri{x) = rjoo 
for ip{p) = ip{p) and the solution ri{x) = Cx for ip{p) = —ip{p) remain valid in the presence 
of impurities. Similar results hold also for intermediate crystal orientations. 

The above analysis shows that in the case of specularly reflecting boundaries the values 
1/(0) and Tjoo are of the same order of magnitude only for particular orientations of the 
normal n within narrow angular intervals A(/)o ~ A^o ~ {io/^{T)Y^'^ around the crystal axes 
Xo, yo, Zq. For other crystal orientations from the equation (|TB]) one has 77(0) ~ Vooio/iiT) <^ 
?7oo, i.e. the d-wave superconducting order parameter turns out to be strongly suppressed in 
the vicinity of the insulating barrier. In the case of diffusive scattering at the interface one 



has g ~ ^0 m^l- Thus in this case the parameter 77(0) ~ ^oVoo/^{T) <^ r/oo for T ^ T^ and 



all crystal orientations. 



At low temperatures the analysis of the proximity effect becomes more difficult. As we 
already discussed for particular crystal orientations the form of the order parameter can be 
described with the aid of symmetry arguments. E.g. for ip{p) = ip{p) the superconducting 
order parameter A(T) does not depend on coordinates whereas for ipiP) = ~'^{p) at any T 
we have A = at the superconductor-insulator interface. For other crystal orientations the 
behavior of the order parameter at T <C Tc can be qualitatively described by the following 
estimate. 

Let us consider the exact equation (|^ and split the frequency range into two intervals: 
\uJm\ -C Aq and {uj^l ^ Ag. The contribution of the first frequency interval to the selfconsis- 
tency equation is small in the parameter cj^/Ao- Therefore for our estimate it is sufficient to 
restrict our consideration to the second frequency interval. For \uJm\ ^ Aq one can neglect 
nonlinear terms in (|^). Then making use of the condition ^o ^ Vx/uJm one gets 



fi{p,x,uj^) = — ■ ■ — + ^ exp(-| \x). (17) 



Taking (0) as an approximate form for /i in the whole frequency interval Um ~ Ag, 
setting ri{x) =const and expanding in powers of the anisotropy parameter /„^>o[V'(p) ~ 
ij{p)]^d^S/ J i)^{p)(fS we find 

This estimate provides correct limits 77(0) = t^oo for ip{p) = ip{p) and r/(0) = for tp{p) = 
—ip{p) and qualitatively describes the low temperature behavior of 77(0) for intermediate 
crystal orientations. It demonstrates that at T <C T^ the value 77(0) is of order rj^o for a 
relatively wide angular interval A0o, A^q ~ 1- The typical length scale at which the value 
ri{x) changes from 77(0) at the boundary to rjoo deep in the superconductor is of order ^q- 

In the case of an ideally transmitting normal metal-superconductor interface D{p) = 1 
quasiclassical propagators are continuous at this interface: g-{x = 0) = g+{x = 0) and 
/_(x = 0) = /+(x = 0). Then imposing the boundary conditions /i(x -^ 00) = f^o, 
A(x -^ 00) = Aoo, /i(— 00) = and assuming that the order parameter is equal to zero in 
the normal metal A(a; < 0) = one can repeat the above analysis and show that at T close 
to Tc the order parameter is again described by the equations (^ and (0), where 

q=( "' f^vp)\, fd^s + ^^^^^ ^^ rmvid'sf X r , .. . . .^2g.-i .^g. 

^ b36C(3)Ty ^ ^^^' "' "^ + 47r3r ! 'il^'^{p)\v,\d^S PJ ^ ^^^ ^ ^ ' ^ ' 

According to this result for NS structures at T ^ Tc the parameter q is of order ^0 for all 
crystal orientations. As before at T — i> the order parameter changes from A(x = 0) to 
A(x = 00) at distances of order ^0 from the NS boundary. To estimate the value 77(0) for this 
temperature interval one can follow the procedure developed in [^ for the case of isotropic 
s-wave superconductors. Then similarly to [^,0 one finds 77(0) ~ 0. 57^00 at T — ^ 0. This 
estimate appears to hold for any crystal orientation. It also agrees with the results of ||3^ 



in which the order parameter of a d-wave superconductor has been calculated numerically 
for a particular crystal orientation. 
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Finally let us note that at T close to Tc and distances x ^ ^o from the interface we also 
expect an additional suppression of the order parameter with respect to that described by 
the equations (|T3|) and (]T5|), (|19|). Indeed, the Ginzburg-Landau equation and its solution 



(i) apply only at distances x :^ ^o from the boundary. For x ^ ^o it is necessary to proceed 
within the framework of a rigorous microscopic analysis |Q. E.g. one can show ||3^ that the 



exact value 77(0) at the boundary between a normal metal and an s-wave superconductors 
is by the factor ^ 1.4 smaller and the exact ratio ?7'(0)/?7(0) is by the factor ^ 1.6 larger 
as compared to the corresponding values which follow from the standard Ginzburg-Landau 
analysis (see e.g. |^). We believe that similar situation takes place also for superconductors 
with anisotropic pairing considered here. 

III. JOSEPHSON CURRENT FOR ANISOTROPIC SUPERCONDUCTORS 

A. Tunnel Junctions 

Let us investigate the dc Josephson effect in tunnel junctions between two d-wave super- 
conductors. Assuming that the junction transparency is small D{p) <^ 1 one can proceed 
perturbatively and expand the boundary conditions (^) in powers of D{p). Keeping only 
the linear terms one gets 

9-{P-) - 9-{P-) = ^I^(P_)(/+(P+)/- (P-) - f-iP-)ft{P+)), (20) 

where the functions f±{p±) and f±{p±) are the anomalous Green functions calculated on 
both sides of the tunnel barrier for D{p) = 0. Substituting this expression into the formula 
for the superconducting current 

JiR) = -2meTj2 f -—^VFiP)gip,R,uJ (21) 

we arrive at the general expression for the Josephson current 

3s = 27reTsin^^ / D{p^)v,{p_)[h4P-)fi+{P+) " f2-{P-)f2+{p+)]^, (22) 

where ip is the phase difference between two superconductors (i.e. the gap is proportional to 
exp[iip)ilj+{p^) and tp-{p_) on the left and on the right respectively), Vx is the Fermi velocity 
projection on the normal to the plane interface. The functions /i, /2 are calculated for real 
V'+(p+), Tp-ipJ). Provided the functions f±{p±,ujm) do not depend on space coordinates 
one can easily evaluate js (|2^) and get 

. ^ . ^ ( A+{p^)A4p_)D{p_)vx{p_) <fS^ 
?« = 27reT sm 09 > / - — -- — . (23) 

tr Jv.>0 , /A2 rr, . ^ + r,;2 . /A2 (f, \ + ,,;2 {2-kYvf 



This result coincides with that obtained in [^. At low temperatures T ^ A it yields 



^,>o |A+(p+)| + |A_(p_)| V|A+(p+)| + |A_(p_)|/ {27TyvF 

(24) 

where K{t) = /g (1 — t^ sin^0)~^/^(i0 is the complete elhptic integral. At A_(p_) ~ 
A+(p^) ~ A(p) and any T we find from ( P5| ) 

This result shows that the Josephson current between identical similarly oriented d-wave 
superconductors at T <C A is proportional to the product A{p)D{p)vx{p) averaged over the 
momentum directions at the Fermi surface. 

With the aid of a standard expression for the normal state resistance of a tunnel junction 

R~N = 2e^ / ^(P->^(P-)7ikr (26) 

Jv^>0 [ZTTJ'^Vf 



one can easily reduce the equations ( p^ and ( pSf ) to the analogous results for conventional 
superconductors provided the momentum dependence of the gap function A is neglected. 
Close to the critical temperature T r^ T^ the expression (l2^) reduces to 



js = — ^ / A4p^)A4p_)Dip_)v.,{p_)-—^. 27) 

As it was demonstrated in the previous section in the presence of a nontransparent 
insulating barrier the order parameter and the Green functions of a d-wave superconductor 
do not depend on coordinates only provided one of the principal axes is perpendicular to 
the barrier plane. For any other crystal orientation the order parameter A as well as g- 
and /-functions vary in space and the results (^)-(^), (0) are no longer valid. To derive 
the corresponding generalization of Eq. (|27|) let us combine the exact formula ( ^21) with the 
results obtained in Section 2. As for a spherical form of the Fermi surface the value D[p) 
depends only on p^, with the aid of Eq. ([131) at T close to Tc we find 



r d?S_ 

js = 87TeTsmipJ2 D{p_y^{p_)I+{p^,uJm)I-{p^,uJm) ,^ .3 . (28) 

^Jv,>o [2'kYvf 

The functions I± read 

1 f°° 2uj 

I±{p±,^m) = -, — i A±{p^, x) exp{-\ — -\x)dx, (29) 

\Vx\ Jo Vx 

X is the distance from the interface. 

The results (|^), (|29|) show that the expression for the Josepson current in d-wave su- 
perconductors is essentially nonlocal in space: the value js is determined by the order pa- 
rameter at distances ~ ^0 from the tunnel barrier. E.g. for a particular crystal orientation 
'4'{p) = ~4'{p) at both sides of the barrier close to Tc we get 

10 



Js = iv'iO)r^^[ ^AP^)^-ipJDipJvlipJ-f^. (30) 

This result demonstrates that even if the order parameter of a d-wave superconductor van- 
ishes at the interface rj{0) = the Josephson current remains finite. The nonlocahty of js 
(p8|), (|29|) results in different temperature dependence of the Josephson current for different 
crystal orientations. E.g. for a particular case ip{p) = ip{p) and T close to T^ from (^) we 
have 

i5Cxr/2(0)(xr,-T, (31) 

whereas for ipiP) = ~4'{p) the result (|30D yields 

js(xr]'\0)cx{T,-Tf. (32) 

To evaluate the Josephson current for arbitrary crystal orientations let us combine the 
results (|^), (jT^) with the equations (^), (^). Then making use of ( |2B| ) and assuming 
D{p) oc pt/p^F (k = 0, 2, ...) H we get 

sin(^ r]l^{k + 2) r f 7C(3), , ,. , 

+ '%i-)'l'-iP-)MP.)"^^^i'S. (33) 



Here g-t are the values of the parameter q (16) on both sides of the interface. The result (^ 



is valid provided the parameters q± are not very large q± -C ^±(^)- According to (|TBp this 
implies that the direction of either one of the crystal principal axes should not be very close 
to that normal to the interface A0o, A^o > {.io/^{T)Y''^. At A0o, A^o ^ (^o/^(^))^/^ the 
expression for js shows a crossover from (^) to (pTf ) . Thus we can conclude that at T close 
to Tc for a wide range of crystal orientations the proximity effect strongly influences the dc 
Josepson current in d-wave superconductors leading to the temperature dependence js (x. 
{j]oo I i{T)Y oc {Tc — Ty. Only if one of the crystal principal axes is (nearly) perpendicular to 
the junction plane this dependence changes and becomes js oc Tc — T. At lower temperatures 
T ^ A(T) the role of the proximity effect becomes less important and the expression ( ^31) 
qualitatively describes the dc Josephson current for a wide range of crystal orientations. 



Sigrist and Rice |T^ suggested the following simple phenomenological expression for the 

^xo ~ Pyo 



Josephson current between two tetragonal superconductors with pi — pi type of pairing 



near Tc. 

Js = Jo sin(^(n^(,^) - n^(j^^,))(n2_(^^^) - n2_(^^^)). (34) 

Here n_|_(j) denotes the projection of the unit vector normal to the boundary plane on the i 
axis. 

Since the basis function ip{p) is not unique for a given pairing symmetry type, the 
dependence of the Josephson current on the orientation of the boundary plane relative 
to the crystal axes cannot be chosen unambiguously from the symmetry arguments only. 
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Eq. (^4]) presents the simplest example for a particular angular dependence of the Josephson 
current consistent with the pairing symmetry. For more complicated cases higher powers of 
^±(x,j/) can appear. 

In the case of diffusive scattering at the boundary the Josephson current does not depend 
on the relative orientation of two superconductors and ri{0) ~ ('Co/'C(^))^^oo [|10l- If ^^ put 
ri{0) ^ , in accordance with (p2|),(p3D at T close to Tc we have jo oc {Tc — T)^. 

As it follows from our analysis the dependence of js on the relative orientation of super- 
conductors becomes important for the specular scattering at the insulating barrier. From 
Eq. (^) one can easily recover the dependence of the Josephson current js on the angle </> 
between the ZQ-axes of two superconductors. E.g. if n±(^x) = 1 and D oc pl/pp we have 

3sR. = ^i^(9 + \cos2^), (35) 

where Ai, A2 are the maximum values of A(p_,_), A(p_). We see that in the case of Eq. 
(P5|) the Josephson critical current is positive for all values of (0-contact). However if, for 
instance, the 2;o-axes for both superconductors are perpendicular to the boundary plane the 
Eq. d^ yields 



js oc cos 20, (36) 

where is now the angle between the xo-axes of the superconductors. In the latter case the 
Josephson critical current turns out to be negative (vr-junction [0) for certain values of 0. 
Thus, rotating one of the crystals around its ZQ-axis one can turn the 0-junction into the 
TT-junction. The latter property remains also for lower temperatures T <^ Tc in which case 
the 0-dependence of the Josephson current is more complicated. 

Note that the dependence ( P^D permits to realize a very simple configuration, for which 
three parts of the same superconductor form a closed circuit with the odd number of n- 
junctions (see Fig. 3). The zo-axes of the grains coincide and are taken to be perpendicular 
to junction planes. The angles between the xo-axes of the first and second as well as the 
second and the third grains are supposed to be equal to 7r/6. Then according to (|36D the 
corresponding junctions are of the 0-type. In contrast the junction between the third and 
the first grains turns out to be of the vr-type because the angle is equal to 7r/3 for this 
junction. 



B. SNS Junctions 

Let us now consider the dc Josephson effect in planar structures superconductor-normal 
metal-superconductor (SNS). In the case of clean s-wave superconductors this effect was 
studied in Refs. [^^. Here we provide the generalization to the case of d-wave supercon- 
ductors with arbitrary orientations relative to the boundaries. 

In order to evaluate the supercurrent in SNS junctions one has to solve the Eilenberger 
equations ([I|) in superconducting and normal regions and match the solutions at NS in- 
terfaces with the aid of the boundary conditions (^). Below we shall consider the case of 
transparent NS boundaries with D{p) = 1 and assume that the order parameter in the 
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normal metal is equal to zero A{—d/2 < x < d/2) = 0, c? is the thickness of a normal layer 
between two d-wave superconductors which occupy two halfspaces x < —d/2 and x > d/2. 
To find the Green functions of superconducting banks we shall use the following ansatz 

f±iP±,x,iu^m) = e*^±(P±) Te^'^±(P±)sgnt;,5(7±(p±,x,cu^), 

fiip±,x,ujm) = e-^^±(P±) ± e-'^^^P±hgnvJg±ip^,x,iuJ, (37) 

g±iP±,x,uj^) = Sg±{p^,x,Um). 

Here (^±(p±) are the phases of the order parameters A±{p_^_). Eq.(^) satisfies the normal- 
ization condition (j^). The ansatz ( |57D is correct for Um <^ ^(p)- The latter inequality in 
turn holds for the parameter region vp/d -C T -C |A-|-(6' = 0)| or T -C vp/d, d ':^ ^q which 
will be considered below. 

Substituting (|37D into (P in the main approximation one obtains 

\v^\VJg±{p,x,Um) = T'2\A±{p^,x)\6g±{p^,x,uJm)- (38) 

The solution of the Eilenberger equations (P in the normal metal is trivial. Combining 
this solution with (P7D and making use of the continuity condition for the Green functions 



at NS interfaces we find 

gN{p,UJm) = 5g^{p^,uJrn) = 6g+{p^,uJm), 

(eiv-(p-) + e'^-'^P-hgnvJg-{p_,uj^))e-^'^^'^/''^ = (39) 

^ gi^+(p+) _ e'^P+^P+hgnvJg+ip^^Um), 

gN is the Eilenberger Green function in the normal metal. Similarly to the case of conven- 



tional superconductors (see e.g. ||33|) the equations (pQ]) yield 



gN{p,uJm) = sgn f ^ tanhj — \ } , (40) 

Z Vx 

where ^p{p) = (p for ip^{pj^)ip_{p_) > and ip{p) = y^ + vr for ip+{p_^_)ip_{p_) < (as before 
the gap is chosen to be proportional to ip-{p_) on the left side and to ex]i{iip)ip+{pj^) on 
the right side of the barrier) . 

Substituting (^0|) into (^If ) we arrive at the final expression for the Josephson current in 
SNS junctions. For vp/d <^T <^ |A±(^ = 0)1 we reproduce the standard result 

is = 6enexp{—2nTd/vF) simp' /md, 



derived before for conventional superconductors p^j38[1 . Here A±{6 = 0) = A{px± = Pf) 
and ip' is the total phase difference between A+(^ = 0) and A_(^ = 0), n = p^p/Sn"^ is 
the electron concentration. The difference between s- and d-wave superconductors becomes 
important in the low temperature limit T <C vp/d. At T —* and d ^ ^o "we get 

J^=4;^(C^iM+C^2b + vr]). (41) 

The function [(p] defines the standard sawtooth behavior oi jsi'p) for s-wave superconductors 
at T = |3§ (see Fig.4) and 



Ci= f cos'^edQ+, C2= f cosHdQ- 
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dQ~^'~ are the solid angle elements on the Fermi sphere for which the functions ip-{'p_) and 
V'+(P+) have equal or opposite signs respectively. The phase dependence of the Josephson 
current in SNS junctions between d-wave superconductors jsi,^) ( PD is presented in Fig. 5a. 
In contrast to the analogous dependence for s-wave superconductors (Fig. 4) it contains an 
additional jump at </? = 0. This jump is due to the presence of an additional phase shift n 
acquired by electrons with momentum directions corresponding to different signs of the gap 
functions in two superconductors. 

We believe that the above mentioned unusual behavior of d-wave SNS junctions can 
be used to provide an experimental test for the symmetry of the order parameter in high 
temperature superconductors. Let us consider a superconducting ring interrupted by an 
SNS junction with the current-phase relation (Pl). Rewriting this relation in the form 
is = Aiip^ A271 respectively ioi < ip < it and — tt < ip < one can easily derive the free 
energy of the system F. In the absence of an external magnetic field we have 

F{I) = ^{P + Kp-2m\I\), (42) 

where / is the current in the ring, L is the ring inductance, /(O) = A2Tt, k = 27cLAi/^q, 
$0 is the flux quantum. This expression is valid for |2L//$o| < l? for larger values of |/| 



the two last terms in ( ^2[ ) are periodically continued with the period ^o/L. The free energy 
of a SQUID with an SNS junction is shown in Fig. 5b for the large inductance limit. It has 
two minima at (p = ±itA2/Ai which correspond to the condition Is{f) = 0. Thus an SNS 
junction between two d-wave superconductors has a twofold degenerate ground state inside 
the interval — vr < (f < tt. This behavior differs from that for tunnel junctions in which case 
the system has only one energy minimum at (y9 = or (y9 = vr. 

Minimizing (|^) with respect to I we find the equilibrium value for the current 

/ = ±/(0)/(l + fi:). 

This result means that an SNS junction described by the current-phase relation (|4l|) always 
induces a spontaneous current in a superconducting ring no matter how small the inductance 
L is. This result differs from that obtained for a ring with a vr-junction |jl2| in which case 
the spontaneous superconducting current can occur only provided L is sufficiently large. 

Without an external magnetic field the ground state of the system is degenerate with 
respect to the direction of the current / fiowing across the ring. This degeneracy is lifted by 
an external magnetic flux $ applied to the ring. In this case the value I in the last two terms 
of (|42|) should be substituted by / + {^/L) and the energies of the two lowest states differ by 
AF = 2J(0)<I'/(1 + k). For k > 1 and Ai ~ A2 we obtain a simple estimate AF ~ <l>o$/I/. 

If one considers a SQUID configuration with two SNS junctions one can easily see that 
the critical current through this system Imax Kiay reach its minimum value not only at 
$/$o = (as in the case of a SQUID with 0-junctions) or at $/$o = 1/2 (as for a SQUID 
with TT-junctions) but at an arbitrary value of $/$o depending on the relation between Ai 
and A2. The dependence Imaxi,^) for a SQUID with identical SNS junctions and Ai > 2A2 
is depicted in Fig. 6. The minimum value of Imax is reached at $/$o = (^1 ~ A2)/2Ai. 
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C. Short Weak Links. 

In addition to tunnel junctions and SNS structures another type of weak hnks between d- 
wave superconductors is of physical interest. Let us consider two superconductors separated 
by an impenetrable insulating barrier with a small orifice of a typical size L '^ ^q. Below we 
shall assume that electrons can freely move (rather than tunnel) through this orifice and put 
its transparency coefficient equal to one D{p) = 1. This model describes various geometries 
(microconstrictions, microbridges etc.) which provide a direct contact between two metals. 
The normal state conductance of such systems depends only on the cross section area of the 
orifice A and is given by the well known expression for the inverse Sharvin resistance 

1/Ro = e^plA/An^. (43) 

In the case of conventional superconductors the dc Josephson effect in this type of weak 
links was studied in details by Kulik and Omel'yanchuk [^. It was found in |2^ that at 
low temperatures the corresonding current-phase relation deviates from the standard sin ip- 
form leading to a somewhat higher Josephson critical current than that for tunnel junctions 
||39|| . Here we briefly discuss the generalization of the theory p8[ for the case of d-wave 
superconductors. 

Following [^] we shall assume that the gap function A is not disturbed in superconduct- 
ing bridges due to the presence of a microconstriction. This assumption is valid everywhere 
except for a narrow region 6r <C C,o close to the oriflce. It is straightforward to check that the 
particular form of A in this region is not important for calculation of the current through 
the oriflce. Therefore without loss of generality (and also for the sake of deflniteness) we 
stick to the same form of the order parameter A(x) in two superconducting bridges as that 
discussed before for the case of tunnel junctions. 

First let us cosider crystal orientations ip±{p) ~ V'±(p) foi' which the value A is (nearly) 
uniform in both superconductors. Then following the procedure P3 one can easily solve the 



Eilenberger equations in superconductors. Matching the Green functions at x = for the 
electron trajectories passing through the oriflce and assuming these functions to be equal 
to the equilibrum ones far from the weak link x -^ ±oo similarly to [Q we obtain the 
expression for the superconducting current through the oriflce Is = Js-A. 

J _o 7-4 V / t;^(p^)A+(p^)A_(p_)siny) <fS^ 

•^ ^ L>o ul + [{ul + ^\{pM< + ^UP-W + A+(p^)A_(p_) cos^ {27^fv, 



m>0 



(44) 



At T > A Eq.dl) reduces to Eq.(0) with D{p) = 1. To analyse the result (^) at lower 
temperatures it is again convenient to introduce the quantity f{p). Then for |A_(_(p_,_)| ^ 
|A_(p_)| similarly to the case of conventional superconductors p8| from (^^ we get 

7. = 2«/ |A_(P-)hmb(p)/2]tauh( l^-'''->l7''^'P)/^l )„.(P-)y^. (45) 



lv:,>o \ 2T J {2tt'^)vf 

Here d'^S' is the element of the Fermi sphere for which the condition |A+(p_(_)| ^ |A_(p_ 



IS 



satisfled. As in |Q the current turns out to be discontinuous at v? = vr. In the opposite limit 
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|A+(p^)| ^ |A_(p_)| and at T = with the logarithmic accuracy we find (0 < (p{p) < 27r) 
I. = 4e£^^.n^(p)|A_(p_)|l„^±(M.ap_)^. for Mp)-.|»,„-^ 

Is = -A7re sgn{ip{p)-7r)\A{p_)\v,{p_) , for |(^(p) - 7r| < In"^ + ' 

(46) 

where (i^S* denotes the element of the Fermi sphere with |A+(p_^)| ^ |A_(p_)|. We see 
that the magnitude of the current jump at ^(p) = tt ( ^G)) is by the factor ~ ln^^(| A+|/| A_|) 
smaller as compared to the case |A+| = |A_| (^). For (p{p) close to vr and arbitrary ratio 
(I A+l/l A_|) the magnitude of the jump reads 

Is = -47re / sgn </. p - ^ , ' . || ' , w. > . P- .. xg " ^7 

Jv,>o \A{p^)\ + \A{p_)\ {2ny^VF 



The integration in ( ^71) runs over the parts of Fermi-surface where (p{p) is close to vr. As 
the function ^p{p) changes its sign on the Fermi-surface an additional jump on the Isiy^) 
dependence takes place at T — > similarly to the case of SNS junctions. 

Let us emphasize again that the result (^) holds only for a homogeneous distribution of 
the order parameter in superconducting banks. Within the same framework an analogous 
result was recently derived by Yip [^. Provided the condition |A(p)| ^ |A(p)| is not 
satisfied the superconducting order parameter depends on the coordinate and the expression 



for Josephson current deviates from (^4[) . In this case after a straightforward calculation 
one can show that at T ^ ^(T) the value Is is given by Eqs. (0), (|29| ) and hence again 
reduces to the result (|33|) with D{p) = 1 {k = 0) and R^ -^ Ro- 

IV. QUASIPARTICLE TUNNELING AND PHASE FLUCTUATIONS 
A. Low Voltage Conductance and I-V Curve 

A possible way to test the symmetry of the superconducting order parameter is to mea- 
sure the I-V curve of a tunnel junction in the limit of low temperature and voltage. In the 
case of isotropic s-wave superconductors at T ^ A only a small number of quasiparticles 
activated above the gap contributes to the junction conductance G. Therefore in the limit 
of small voltages we have G oc exp(— A/T). At T = no quasiparticles exist above the gap 
and the current across the junction is equal to zero / = provided the externally applied 
voltage V does not exceed the value A/e for NS junctions and 2A/e for SS junctions. Below 
we shall show that in the case of d-wave symmetry of the order parameter the I-V curve of a 
tunnel junction is entirely different in the corresponding temperature and voltage intervals. 

Let us assume that the time-independent external voltage V is applied to the tunnel 
junction between two metals. Then expressing the current in terms of the Green function 
of the system and making use of the boundary conditions (^) in the lowest order in D after 
a standard calculation (see e.g. [^) one easily finds 
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•^ ^^JUv ^^^^(2^) ~^^'^^( ^ 2T ) 3+i^ ~ ^^^P+'>3'-(^^P-)^^)^^(P^)^(P 



£S. 



(48) 



Here jat is a dissipative contribution to the current across the junction and g'j^{e,p^) are the 
normahzed densities of states of two metals in the vicinity of a tunnel barrier. In the case 
of d-wave superconductors for A(p) = A(p) we have 



2 



^Up±) 



Let us first calculate the I-V curve of an NS junction. Setting A_ = and A+ = A(p) 
and substituting (E^) into (^8]) we obtain at T = 



2e / [{eVf - ^\p)Yl'Q{eV - |A(p)|)D(p)^;.(p)-^. (50) 



jN = 2e I [{evy - /\-{p)r-<d{eV - \/\{p)\)V{p)v^{p)-—^ 
Jv^>o [ZiryvF 

The equation (^Dj) defines the dissipative current across the tunnel junction for crystal ori- 
entations with A(p) ^ A(p). For other crystal orientations the superconducting density of 
states in the vicinity of a tunnel barrier deviates from (|^). Nevertheless - as in the case 
of a dc Josephson current - at T = the result (|50D remains to be valid apart from an 
unimportant numerical factor of order one. Below we shall neglect this factor and apply 
the result ([50| ) to any crystal orientation. Then choosing the order parameter in the form 
A(p) = Ao(p^Q - pU from Eq. (0) we have 

j^ = ^ f" d(j) r de sine coseD(e){(eVf - Allsin^ e cos^ (j) 
Att-^ Jo Jo 

- (sin e sin (f) cos ^o - cos 6 sin ^o)T}^^^- (51) 

Here 60 is the angle between the vector n normal to the junction plane and the crystal axis 
Zo, Aq is the maximal value of A(p). 

It is easy to see that - in contrast to the case of s-wave superconductors - the current 
(pT|) does not vanish even for eV ^ Aq. In the latter limit the main contribution to Jn 
comes from quasiparticles with the momentum directions close to the directions for which 
the order parameter A(p) is equal to zero. The integral over these momentum directions 
can be in turn splitted into two terms 

JN = jm+JN2- (52) 

The first term j^i is defined by the integral over the values 6 close to 60 and all values of (p 
from to 27r or, in other words, over the momentum values Pz^ ~ Pf- With the logarithmic 
accuracy the corresponding integration in (^) yields 



ev,ieo)Dieo)ieV)yj, ^ A 







jNi = ^^ — T In --. (53) 
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The second term Jn2 comes from the integration over momentum directions close to the 
hnes PxQ = ±Pj/o- In the vicinity of these hnes the gap function is A(p) = {p2/pF)^oh{pi) 
where pi is the coordinate along the line of zeros and p2 the one in perpendicular direction. 
In our case h{pi) = 2sm9', where 9' is the angle between p and Zq. Integrating over p2 we 
obtain 



JN2 



eVfpF 
Stt^AqVf 



I ^^irr^dp]^^^. (54) 



Comparing the results (|53|) and (|5^) one can conclude that in the limit eV <^ Aq the current 
JNi dominates for crystal orientations 7r/2 — 16*01 S> ln~^''''^^"'^^(Ao/eV^). E.g. for 6*0 = and 
D{e) occos^^ {k = 2) we get 

For 7r/2 — |6'o| <C ln~ ' (Ao/e^) the axis zq nearly coincides with the junction plane and 
the term j^i becomes small. In this case the current j^ is given by the term jn2 (|5^)- For 
0Q = ±7r/2 and eV < Aq it yields 

The zero temperature I-V curves for an NS tunnel junction with D oc p'^/p'j? are pre- 
sented in Fig. 7 for two particular crystal orientations (one of the principal axes xq or zq is 
perpendicular to the barrier plane). At low voltages eV <^ Aq these curves follow the results 
(^), (0) (improving the logarithmic accuracy of (pS)): ln(Ao/eV) -^ ln(2.4Ao/eV)). At 
higher voltages eV ~ Aq the I-V curves shows a smooth crossover to the standard Ohmic 
behavior. 

The I-V curve of a tunnel junction between two d-wave superconductors can be calculated 
analogously. Substituting (^) into (^) for the case of identical superconductors at T = 
we find 

r / .ey-|A(p+)| uj{eV-uj)duj \ ,„ , ,. , d^S. 



v.>o\J\Aip_)\ (u;2-A2(p_))i/2((eV-^)2_A2(p^))i/2; ^^-^ "■^^-^(27r)3y^- 

(57) 

The integration in ( |F7| ) is made over that parts of Fermi-surface where eV — |A+(p^)| > 
|A_(p_)|. In a general case the zero lines of the order parameters in two superconductors 
do not coincide. Assuming that the angle x between these lines in the point of their inter- 
section Pj obeys the condition e\^/Ao|/i(Pj)|x ^ l? we can easily obtain the leading order 
contribution to the quasiparticle current for eV ^ Aq: 



eVx{Pi)D{Pi)pUeVY 
2ATTVF\h+{pi+)h4pi_)\ I sin x| Ag ' 






In order to find the total current it is necessary to sum up the contributions from all inter- 
section points. Then one obtains 



e2y3 
Jn = a-^—r^- (59) 



Here the factor a keeps track on the particular relative orientation of superconductors and is 
of order one for most of such orientations. This factor vanishes only provided Vx{pi)D{pi) = 
0, i.e. if the intersection points coincide with the poles of the Fermi surface and the Zq axes 
of both superconductors are in the junction plane. 

If both superconductors are oriented identically we again arrive with the aid of ( |57| ) at 
the expression for the current j^ defined by the expressions (^31)- (^61) multiplied by the 
numerical prefactor 



71 Jo \1 — UJ J 

where as before K{t) is the complete elliptic integral. In this case the zero lines of A_(p_), 
A_|_(p_,_) coincide if they are drawn on the same Fermi- surface. 

The zero temperature I-V curves for the junction between two d-wave superconductors 



calculated numerically from the equation (|57|) are presented in Fig. 8. The upper and the 



lower curves were calculated assuming that respectively zq and xq axes of both supercon- 
ductors are perpendicular to the boundary plane. Orientation of Xq and yo axes of two 
superconductors is identical for the upper curve whereas for the lower curve their yo axes 
constitute the angle it/2 between each other. Again in the low voltage limit the numerical 
curves agree well with the analytic results ([59|), (^) and allow to define the corresponding 
numerical prefactors. E.g. the factor a in Eq. ( p9D is found to be equal to a ~ 0.35 for the 
above crystal orientation and the logarithmic accuracy of Eq. (|55D can be improved by a 
substitution ln(Ao/e1/) -^ ln(3.9Ao/e1/). 

For larger voltages of order Aq the differential conductance G = dl/dV has a maximum 
which exact position depends on the relative crystal orientation. E.g. for the two curves 
presented in Fig. 8 the conductance G{V) reaches its maximum respectively at eV ~ I.OSAq 
and eV ~ I.QTAq. 

In order to understand the difference in the maximum positions by nearly a factor 2 
for these two configurations let us consider the case A_|_(p_,_) ^ A_(p_) and carry out the 
frequency integration in Eq. (^Tj) . Then we obtain 



dl _ J^^^,k-'[E{k){l + P) - (1 - e)K{k)]D{p_)vxip_)d^S. 

dv'''' 2J^^^,D{p_)vx{p_)d^S^ ' ^""^ 

where fc(p) = {{eVf - {2 A{p)fy/^/ eV and E{k) = /o''/^(l-A;2 sin^ 0)1/2^/0. The momentum 
integration in (^) is carried out over that parts of Fermi surface where eV > 2|A(p)|. It is 
easy to see that the maximum of dl/dV takes place at a certain effective gap value which 
(due to the presence of the factor D{p_)vx{p_) in the integrand (|60|) ) is mainly defined by 
the value of gap function in the direction normal to the junction plane. Since for the lower 
curve of Fig. 8 for both superconductors the gap function is equal to Aq in this direction 
the maximum of dl/dV takes place at eV ~ 2Ao. For the upper curve the gap vanishes 
along the direction normal to the interface. Accordingly the maximum has a much smaller 
amplitude and takes place at lower voltages eV ~ Aq. 

Note that the behavior of the low voltage conductance G oc V^ has been detected in 



recent experiments with SS tunnel junctions |T9|. This behavior is in a good agreement with 
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our theoretical predictions (pHD and (|59[). Also for higher voltages our results qualitatively 
agree with those reported in [jl9|. 

The dependence G ocV^ for the low voltage conductance of a tunnel junction between d- 



wave superconductors has been also discussed in a recent paper by Won and Maki [^ within 
a different theoretical framework. They evaluated the quasiparticle current by means of 
the standard tunneling Hamiltonian approach assuming that tunneling matrix elements are 
independent of the momenta of tunneling electrons and making use of the expressions for the 
superconducting densities of states averaged over all momentum directions. This approach 
yields the results which are independent of relative orientation of two superconductors. 
Although it appears to be quite difficult to justify such an approach microscopically we 
believe that it might work - at least qualitatively - for diffusive SS boundaries. However 
it clearly fails for specularly reflecting boundaries in which case the quasiparticle current 
essentially depends on the relative orientation of d-wave superconductors. 

Very recently the case of specularly reflecting boundaries has been independently studied 
by Bruder, van Otterlo and Zimanyi [^. These authors also proceeded within the tunnel- 
ing Hamiltonian approach completed by a phenomenological assumption about the angular 
dependence of the tunneling matrix elements. For identically oriented superconductors with 
zq axes being in the barrier plane they also arrived at the result j^ oc V^ which agrees 
with our results (^^, (|56D. However for misoriented superconductors at T = a vanishing 
subgap current I{eV <^ Aq) = has been found in |^. In contrast our results (0), (|54D 
and (|58[) demonstrate that even at T = the subgap current does not vanish for all crystal 



orientations with Vx{pi)D{pi) ^ |4^, Pi is value of the electron momentum at the inter- 



section point of the nodal lines for the two order parameters A+ and A_. The origin of 



this disagreement lies in the fact that the authors |4^ considered the case of a long cylin- 
drical Fermi surface whereas the analysis developed here is based on a picture of a (nearly) 
spherical Fermi surface. As under certain restrictions the Eilenberger formalism can be also 
applied to superconductors with nonspherical Fermi surfaces [^] the corresponding gener- 
alization of our approach can be easily provided. E.g. Eqs.(^), (|58|) remain valid for the 
metals described by the dispersion law e{p) = {pl^ + Py^)/2mi + pl^/2m2- In order to find 
the expressions for h^ and x iii ( |58D one should rewrite the function A_|_(p_,_) in terms of the 
variable p_ and draw the nodal line of this function on the Fermi-surface of the "-"-metal. 
Then it is easy to see that the prefactor a in Eq.(|5^) turns out to be of order a ~ (mi/m2)^^^. 
It appears that the results of |^ correspond to the limiting case mi/m2 -^ for which the 



prefactor a in (^) (and thus the current Jn) vanishes at T = and small V. Physically this 
situation corresponds to ideally two-dimensional character of the electron motion in xq — yo 
planes whereas in the 2;o-direction this motion is totally suppressed. If one allows for jumps 
of electrons between such planes the ratio mi/m2 differs from zero and the result Jn oc V^ 
remains valid. It is easy to check that the validity condition for this result x ^ e\^/Ao 
remains unchanged for any nonzero mi/m2 whereas the result Jn ~ V^ holds only in the 
narrow region of misorientations of two superconductors x "^ (e^/Ao)(mi/'m2)^''^. 

AT low but finite temperatures T ^ Aq the results obtained above for the case T = 
remain valid for not very small voltages eV ^ T. In the opposite limit eV <C T the main 
contribution to the current Jn comes from quasiparticles thermally activated above the gap. 
In the limit eV -C T -C Aq for the NS junction we find 
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jN = G{T)V, (61) 

where the hnear conductance G{T) of a tunnel junction between misoriented superconduc- 
tors is 

G{T) ~ TyRjvAl. (62) 

Analogously for the crystal orientations described by the equations (|53|) and (0) one gets 
respectively G{T) ~ {T/RnAq) ln(Ao/T) and G(T) ~ {T/RnAq). 



B. Effective Action 

Finally let us briefly demonstrate how the above results can be generalized to take into 
account thermodynamic and quantum fluctuations of the phase difference ip across the tunnel 
junction between d-wave superconductors. The grand partition function of this junction can 
be expressed in terms of the path integral over the (^-variable (see e.g. 



Z ~ Jv^{T)exp{-S,ff[ifir)]), (63) 

r is the imaginary time variable which changes from to /5 = 1/T. To evaluate the effective 
action functional Seff[(p] we make use of the approach developed in Refs. P6|J^ which 
allows to recover S^jf [ip] from the expression for the kernel of the current density operator 
j{r, t) by means of the integration over the effective "coupling constant" A. In our case the 
corresponding formula reads 

S[ip] = f'dX [^ dMT)j[Xip{T)]A/2e, (64) 

JO Jo 

where j[(y9(r)] represents the current density through the junction. In the interesting limit 
of low frequencies the expression for the supercurrent js reduces to the standard Josephson 
relation 

js = jo sin ip{T), (65) 



whereas the kernel for the quasiparticle current operator has the form |^ 



j^[y.(r)] = 2e|^''rfr'a(r-r')sin(^^^^^^'). (66) 



Combining (M)-(|66D and also taking into account the charging energy term one immediately 



arrives at the AES effective action |^| 

'"' = i!^^ii(B' - ^■'"^»''^)i - T'^^r*'"'^ - o "< '''"' 7'''' )- («^) 

C is the junction capacitance and Ej = JQA/2e is the Josephson coupling energy which can 
be positive or negative depending on the relative crystal orientation. The particular form of 
(y{T) depends on the form of the I-V curve in the limit of small V. E.g. making use of 
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it is easy to show that for I oc V^ one has a(r) oc r '^. More precisely, combining (|59|) with 
(|66|) we obtain at T — i> 

a(r) = Sa/Tre^i^TvAgr^ (68) 

According to the above analysis this result holds for most of crystal orientations. For 
the orientations describrd by the I-V curves (|55D and (|56D we find respectively a(r) ~ 
ln(AoT)/e^-RArAoT'^ and a(r) ~ l/e^Rj\fAQT^. The latter dependence has been also obtained 

mil. 

We believe that the above results might be helpful for a quantitative description of ther- 
modynamic and quantum properties of Josephson junctions and granular arrays composed 
by d-wave superconductors. 

V. DISCUSSION 

The microscopic analysis of the charge transport in tunnel junctions and weak links 
formed by d-wave superconductors allows to encover several interesting features of such 
systems. We demonstrated that the order parameter of a d-wave superconductor can be 
essentially suppressed in the vicinity of the insulating boundary depending on its orientation 
relative to the principal crystal axes of such a superconductor. This proximity effect can 
in turn strongly influence the Josephspn current between two superconductors and becomes 
particularly important at T close to Tc. In the latter case the temperature dependence of the 
Josephson critical current jo varies from js (xTc — T for a homogeneous order parameter in 
superconducting bulks (i.e. if one of the principal crystal axes is nearly perpendicular to the 
junction plane) to jo oc (Tc — T)'^ for other crystal orientations. The results of our calculation 
show a significant dependence of the Josephson current on the relative orientation of the 
superconductors and are consistent with those of recent experiments |T^-|T3 which indicate 
the possibility of d-wave pairing symmetry in HTSC compounds. 

The current-phase relation for SNS junctions and short superconducting weak links es- 
sentially deviates from the standard Josephson relation js = jo sin (p in the low temperature 
limit. In the case of d-wave symmetry of the order parameter at T = the current-phase 
relation js{f) for SNS junctions shows an additional jump (as compared to the case of 
s-wave pairing) at the point (f = 0. Accordingly the superconducting coupling energy for 
such junctions E{(f) (in contrast to tunnel junctions) has two degenerate minima within 
the phase interval — vr < (f < tt (see Fig. 4b) which correspond to two different stable zero 
current states. Positions of these minima do not coincide with (f = ot (p = it (as for tunnel 
junctions) but can be located at any point inside the interval — vr < ip < tt. Thus in the 
analogy with 0- and vr-junctions one can say that the systems in question provide an inter- 
esting example of "whatever-junction". Being included into a SQUID ring such a junction 
induces a spontaneous superconducting current no matter how small the ring inductance is 
and yields to further features different from those for SQUIDs with tunnel junctions. We 
believe that these features could be used as an additional test for the pairing symmetry in 
HTSC. 

An additional information about the form of the superconducting order parameter and 
the density of states is contained in the expression for the quasiparticle tunneling current. 

22 



We evaluated this current for tunnel junctions between a normal metal and a d-wave super- 
conductor as well as between two d-wave superconductors in the low temperature limit. The 
corresponding I-V curves show zero-bias anomalies of the type Jn oc V"^, Jn oc V'^\ia{l/V) 
or j]\[ oc V^ depending on the junction type and relative crystal orientations. The latter 
dependence agrees well with the experimental results |19|. At larger voltages the differential 
conductance of SS junctions has a peak (also detected experimentally |l^,|20|) which position 
also depends on the relative crystal orientation. 

We would like to point out that one can in principle provide an example of a vr-junction 
not only between d-wave superconductors but also between s-wave superconductors with 
multisheet Fermi surfaces and different signs of the order parameter on different sheets [^ 



From this point of view an experimental confirmation of vr-j unction-like properties of HTSC 
compounds yet cannot completely exclude s-wave pairing. On the other hand, for special 
types of Fermi surfaces the anisotropic s-wave order parameter can be equal to zero for 
certain momentum directions not only due to the symmetry reasons (as it would be in the 
d-wave case). Therefore the low temperature measurements of a quasiparticle tunneling 
current rather can be considered as a "gaplessness" test than really distinguish between s- 
and d-wave types of pairing. Bearing all that in mind one can conclude that it is quite 
important to combine both dc Josephson effect and quasiparticle tunneling measurements 
for the same tunnel junctions. Although even demonstration of combined vr-j unction-like 
and gapless properties of such systems formally cannot yet exclude other than d-wave types 
of pairing it would strongly favour the possibility of d-wave pairing in HTSC. 

The results derived here are not specific for HTSC compounds and can be also apphed 
to other types of unconventional superconductors, like heavy fermion superconductors. Our 
analysis holds for an arbitrary form of the Fermi surface. For the sake of definiteness some 
limiting results were derived for the case of a spherical Fermi surface. The latter is by no 
means restrictive for any of the concusions reached in the present paper. E.g. our results 
also remain valid for the case of a (nearly) cylindrical Fermi surface which appears to be 
more relevant for several HTSC compounds. The modification of our results for the latter 
case reduces to an effective renormalization of the junction normal state resistance. 

It is important to emphasize that our analysis is completely based on the microscopic 
theory and does not involve model assumptions which are inevitably present in the tunneling 
Hamiltonian approach. Within the latter approach the correct dependence of the current on 
the momentum directions and crystal orientations (important for d-wave superconductors) 
usually cannot be recovered in a unique way and the validity of the final results rather 
depends on physical intuition of the authors than it is controlled by the method itself. 

After this work has been already completed we became aware of the paper [^] where the 
Josephson current between superconductors with mixed s-|-d symmetry of the order param- 
eter has been analysed within the tunneling Hamiltonian approach. The results obtained in 
this paper differ from ours. The source of this discrepancy lies in an incorrect phenomeno- 
logical expression for the superconducting current across the Josephson junction used in 
^ . This expression does not account for a proper angular dependence of tunneling matrix 



elements and therefore leads to irrelevant results. 
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FIGURES 

FIG. 1. The relative orientation of the principal crystal axes xq, i/q, zq and the vector n normal 
to the boundary plane. 

FIG. 2. The parameter q (in units of IGi'^vp/Tc) as a function of crystal orientation relative to 
the boundary plane. 

FIG. 3. An example of the circuit which contains the odd number of tt -junctions. 
FIG. 4. The sawtooth function [p] of Eq.(41). 
FIG. 5. The current-phase dependence (41) (a) and the energy (b) of an SNS junction between 

d-wave superconductors at T = 0. 

FIG. 6. The maximum current Imax through the SQUID with two SNS junctions (described 

by the current-phase dependence (41)) as a function of the external magnetic flux $. 

FIG. 7. The I-V curves for a tunnel junction between a normal metal and a d-wave supercon- 
ductor at T = 0. The results are for the crystal orientations the axis zq is perpendicular to the 
junction plane (upper solid curve) and coincides with this plane (lower solid curve). Ohmic I-V 

curve is shown by a dashed line. 

FIG. 8. The same curves for a tunnel junction between two d-wave superconductors. The upper 

solid curve corresponds to similarly oriented superconductors with their axes zq being perpendicular 

to the junction plane. For the lower solid curve the xq axes of both superconductors were taken to 

be perpendicular to the junction plane whereas their t/q axes were rotated by the angle 7r/2 with 

respect to each other. 
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